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We study the electric conductivity of the vacuum of quenched SU (2) lattice gauge theory induced 
by the magnetic field B as a function of the bare quark mass m. The conductivity grows as the 
quark mass decreases. Simplest power- like fit indicates that the conductivity behaves as B/s/m. 
We discuss the implications of this result for dilepton angular distributions in heavy ion collisions. 
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Heavy ion experiments at RHIC have found an evi- 
dence [1] for the so-called Chiral Magnetic Effect (CME) 
[2] in quark-gluon plasma. The essence of the effect is the 
generation of electric current along the direction of the 
external magnetic field in the background of topologically 
nontrivial gauge field configurations. Experimentally, the 
effect manifests itself as the dynamical enhancement of 
fluctuations of the numbers of charged hadrons emitted 
above and below the reaction plane in off-central heavy- 
ion collisions. In noncentral heavy-ion collisions, very 
strong magnetic field acting on hadronic matter is cre- 
ated due to the relative motion of the ions themselves. 

Recently, this enhancement of fluctuations was also 
shown to be related to the electric conductivity of the 
vacuum of the gauge theory [3]. Namely, it was found 
that an external magnetic field applied to the confin- 
ing vacuum of quenched non-Abelian lattice gauge the- 
ory induces nonzero conductivity in the direction of the 
field. In other directions, the vacuum remains an insu- 
lator. This phenomenon was called "electric rupture by 
a magnetic field" , by analogy with similar phenomena 
known in condensed-matter physics [4] . Electric conduc- 
tivity in the presence of an external magnetic field was 
also calculated in holographic models [5, 6]. 

However, it turned out that theoretical and experimen- 
tal relations between the CME and the conductivity in- 
duced by the magnetic field are far from trivial [7, 8] . In a 
recent paper [7] , it was suggested that massless excitation 
which is responsible for both the CME and the nonzero 
conductivity is the "chiral spiral" , however, no qualita- 
tive predictions for the conductivity were obtained. In 
experiment, the induced conductivity is responsible for 
angular distributions of soft leptons rather than hadrons 
[3] (see also [9]). 

Here we perform the simplest possible check of the 
"chiral" nature of the conductivity induced by the mag- 
netic field: we study its dependence on the bare quark 
mass. As the nonzero quark mass leads to transitions 
between left- and right- handed states, if the imbalance 
of chirality is a necessary condition for the induced con- 
ductivity, it should decrease with the quark mass. Such 



decrease can, for example, serve to distinguish between 
ordinary and strange hadron matter produced in heavy- 
ion experiments [10]. 

The technical setup of our simulations is the same as 
in [3]. Namely, we extract electric conductivity from 
the Euclidean correlator of two vector currents ji (x) = 
q(x)jiq (x): 

Gh(t)= J d 3 x(.n (0)j j (x,r)) = 

= J -^K(w,T)pij H , (1) 
o 

where p^ (w) is the corresponding spectral function, 
K(w,t) = % ^j^f^ [11, 12], T = (Nta)- 1 is the 

temperature, a is the lattice spacing and N t is the size 
of the lattice in the time direction. The phenomenon of 
induced conductivity was found only in the confinement 
phase [3] , therefore in the present paper we also limit our 
studies to the case of small temperatures. 

The underdefined linear system (1) is solved using the 
Maximal Entropy Method [12, 13]. The Kubo formula for 
the electric conductivity yields the conductivity in terms 
of the spectral function in the limit of zero frequency 
[12, 14]: 
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The correlator (ji (x) jj (y) ) is measured on the equi- 
librium ensemble of quenched SU (2) gauge fields with 
lattice tadpole-improved action: 

( Ji (x) 3 'j (v)) = {Q (x) liq (x) q (y) jjq (y) ) = 

-s 
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where is the non-Abelian gauge field with the ac- 
tion Sym [Afj] and V \A^] is the lattice Dirac operator. 
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For V [A^] we use Neuberger's chirally invariant overlap 
Dirac operator [15]. We implement the Shifted Unitary 
Minimal Residue method (SHUMR) described in [16] to 
find the quark propagator. For each value of the quark 
mass and the magnetic field we have used 30 configu- 
rations of gauge fields on 14 4 lattice with lattice spacing 
a = 0.102 fm. In order to study the finite- volume effects, 
we have also repeated some calculations for 50 configura- 
tions on the 20 3 x 14 lattice with the same lattice spacing. 
Uniform magnetic field is added to the Dirac operator by 
substituting su (2)-valued vector potential with u (2)- 
valued one A ' ui j — A^ij + 1/2 Fp V x v Sij. The magnetic 
field is directed along the z axis. In order to account for 
periodic boundary conditions we introduce an additional 
boundary twist for fermions [17, 18]. 

The correlators Gij (r) for various components of the 
currents and for different magnetic fields and quark 
masses are plotted on Fig. 1. As in the case of small 
quark mass [3] , at larger quark masses the correlators of 
the components of the current parallel to the magnetic 
field decay much slower than for the transverse compo- 
nents. As the bare quark mass m q increases, the cor- 
relators also decay faster, however, relative asymmetry 
between longitudinal and transverse components remains 
significant. 

A remark is in order concerning the applicability of 
the expressions (1) and (2) for the lattices which have 
equal temporal and spatial extend. While such lattices 
correspond to zero temperature in the standard lattice 
lore, one can still apply the expressions (1), (2) with a 
small but finite temperature T = (Nta) 1 = 140 MeV = 
0.45 T c , where a is the lattice spacing. This, however, 
might lead to quite large finite- volume effects. In order 
to show that finite volume does not affect significantly 
the vector- vector correlator (1), on Fig. 2 we compare 
the correlators (1) calculated on the 14 3 x 14 and 20 3 x 14 
lattices at magnetic field strength \JqB = 0.45 GeV. One 
can see that the correlators agree within error bars for the 
currents both parallel and perpendicular to the magnetic 
field. Therefore one can expect that our numerical results 
are the same as those obtained on the 14 x L 3 lattices 
with L 14. Thus for the correlators of vector currents 
finite-volume effects are smaller than the statistical er- 
rors in our simulations, and there is no need to use large 
anysotropic lattices. 

Now we analyze the conductivity extracted form these 
correlators using the expressions (1), (2). The depen- 
dence of the conductivity in the direction of the magnetic 
field on field strength and quark mass is illustrated on 
Figs. 3 and 4. At zero magnetic field the conductivity is 
equal to zero within error range for all quark masses. As 
the magnetic field is turned on, the conductivity in the 
direction of the field grows, but remains equal to zero 
for transverse directions. For larger quark masses the 
conductivity in the direction of the magnetic field slowly 
decreases. Transverse components of the conductivity 



tensor remain equal to zero. 

In order to describe our data by some phenomcnolog- 
ical expression, we have fitted the conductivity in the 
direction of the field as a function of both field strength 
and the quark mass by the simplest power-like law: 



(m q ,qB)~mZ a (\qB\y 



(4) 



The fit yields the following values of the parameters a 
and /3: 



a = (0.45 ± 0.06) 
P = (1.1 ± 0.2) 



(5) 



The relative errors of the values of a, f3 are within 20%. 
We can conclude, therefore, that with a good accuracy 
the dependence of the induced conductivity on the mag- 
netic field and the quark mass can be described by the 
following simple expression: 



\lB\ 



(6) 



with some constant A. Such fits are shown on Fig. 4 by 
solid lines. It is interesting that the phenomenological 
expression (6) is non-analytic with respect to the field 
strength and the quark mass. Moreover, if the conduc- 
tivity is indeed given by (6), it diverges at zero quark 
mass. 

Finally, let us relate the conductivity which we have 
measured to experimental observables. Dilepton emission 
rate from either cold or hot hadronic matter is given by 
a general expression [19, 20]: 
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where p\ and P2 are the momenta of the leptons, q = 
P1+P2, L^ u = ((pi -p 2 + m 2 ) rf v ~p^p v 2 -p%p\) is the 
dilepton tensor (rj tJ,L/ is the Minkowski metric), m is the 
lepton mass. If the electric conductivity is nonzero in the 
direction of the magnetic field, for sufficiently small p\, p 2 
one has Pij (q) a^q/T ~ BiBjq/ (\B\T). When the 
mass of the dilepton can be neglected and we are in the 
rest frame of the dilepton pair, one has p\ = —p2 = pn, 
q = ^2p, Oj and the spatial components of the dilepton 

tensor are L ZJ — p 2 (Sij — n, rij ) . The emission rate is 
therefore proportional to 
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(0), (8) 



where 9 is the angle between the spatial momentum of 
the outgoing leptons and the magnetic field. Therefore, 
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FIG. 1: Current-current correlators at different quark masses. Above on the left: correlators of the longitudinal components of 
the current at \/qB = 0.45 GeV and different quark masses. Above on the right: correlators of the longitudinal and transverse 
components of the current at ■s/qB = 0.45 GeV and different quark masses. Below: correlators of the longitudinal and transverse 
components of the current at m q = 100 MeV (on the left) and at m q = 540 MeV (on the right) and different magnetic fields. 
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FIG. 2: A comparison of current-current correlators (1) cal- 
culated at 14 3 x 14 and 20 3 x 14 lattices at equal values of the 
lattice spacing and magnetic field strength \/qB — 0.45 GeV. 
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FIG. 4: Induced conductivity in the direction of the magnetic 
field as a function of the bare quark mass at different magnetic 
fields. 
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FIG. 3: Induced conductivity in the direction of the magnetic 
field as a function of field strength at different quark masses. 

there should be more soft dileptons emitted perpendic- 
ular to the magnetic field than parallel to it. Accord- 
ing to (8), the multiplicity of soft leptons should grow 
linearly with the induced magnetic field. This result, 
however, should be taken with care, since our data are 
obtained for quenched theory at sufficiently low tem- 
peratures, below the confinement-deconfinemcnt phase 
transition. In practice, of course, one could also expect 
numerous soft dileptons coming from the thermalized 
quark-gluon plasma formed in the collision. Neverthe- 
less, "anisotropic" soft dileptons described by (8) could 
be in principle found in heavy-ion experiments such as 
RHIC, FAIR or NICA. The mass dependence (6) of the 
anisotropic part (8) of the lepton production rate could 
be used, for example, to distinguish between ordinary 
and strange nuclear matter formed during the heavy- 
ion collisions. Strangeness enhancement in the vicinity 
of the confinement-deconfinemcnt phase transition [10] 
might be a relevant situation. 
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